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Abstract 

Quasinormal modes provide valuable information about the structure of space- 
time outside a black hole. There is also a conjectured relationship between the 
highly damped quasinormal modes and the semi-classical spectrum of the horizon 
area/entropy. In this paper, we show that for spacetimes characterized by more 
than one scale, the "infinitely damped" modes in principle probe the structure of 
spacetime outside the horizon at the shortest length scales. We demonstrate this 
with the calculation of the highly damped quasinormal modes of the non-singular, 
single horizon, quantum corrected black hole derived in 



1 Introduction 



During the last decade, a great deal of attention has been given to the "ringing" 
or quasinormal mode (QNM) spectra of black holes. QNM's are the damped fun- 
damental frequencies at which black hole spacetimes resonate when perturbed. 
The long-lived QNM's are in principle observable in the gravitational waves emit- 
ted during the ring down of a black hole formation. The highly damped modes, 
while not in practice observable, have in recent years become a subject of interest 
due in part to a conjecture by Hod[T] relating them to the spacing between semi- 
classical area/entropy eigenvalues of the quantum black hole. Hod's argument 
was based upon earlier numerical results for the highly damped QNM's of the 
Schwarzschild black hole and the Bohr Correspondence Principle. 

Our purpose here is not to settle the controversy over this intriguing conjec- 
ture. Whether or not one subscribes to a specific relationship between highly 
damped QNM's and the black hole area spectrum, the impressive body of calcu- 
lations in the literature makes it clear that there is a deep connection between 
the black hole QNM's and the small scale structure of the black hole spacetime. 
The nature of this connection was at first considered a bit mysterious due to the 
fact that the calculation of the asymptotic QNM's required analytic continuation 
of the exterior solution to the entire complex r-plane. This generated questions 
along the lines of: "Why would the vibrational modes exterior to the horizon care 
about the singularity at the origin?" An answer to this was already implied in the 
paper of Andersson and Howls[2] in the context of the highly damped QNM's of 
Reissner-Nordstrom (R-N) black holes. It had been observed that in the highly 
damped limit, the real part of the QNM's of the R-N black hole approached ln(5) 
when the charge q was taken to zero. This appeared to be in contradiction to 
the fact that for the Schwarzschild solution (which of course is the q — > limit 
of the R-N solution) one obtains ln(3). Andersson and Howls suggested that the 
explanation for this apparent contradiction lay in the existence of two separate 
scales in the R-N problem, both in addition to the scale set by the ADM mass 
of the solution. They argued that care had to be taken as to the order in which 
the various scales were taken to zero or infinity. Specifically, the q — > limit 
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does not commute with the \ui\ — > oo limit, where Ui is the imaginary part of the 
complex QNM frequency oj which is also the damping rate. If the former limit 
is taken first, one obtains ln(3) whereas if the latter limit is taken first, ln(5) is 
the result. They noted therefore that while ln(5) is the correct R-N result for 
very high damping, there should be an intermediate range of damping for which 
the Schwarzschild value of the real frequency is correct. Order of magnitude 
arguments suggest this range in four spacetime dimensions (4 — d) to be: 

1 « GM\u\ « (GM) 4 / q A , (1.1) 

where M is the black hole mass and G is Newton's gravitational constant. Note 
that \u\ ~ in the large damping limit. This observation was verified explicitly 
in [3J. A similar situation also appears in Kerr black holes with small angular 
momentum. In the case of Kerr spacetime, it was observed by Keshet and Hod[l] 
that in the highly damped limit, the real part of the QNM frequency approached 
zero when the angular momentum per unit mass a was taken to zero. This is 
again in contradiction to the fact that for the Schwarzschild solution, which is the 
a —7- limit of the Kerr solution, one obtains ln(3). This issue was resolved by 
Green, Mulligan and one of us in [5], where it was shown that in the intermediate 
range of 

1 « GMU « — — (1.2) 

a 

the Schwarzschild value of the real part of the QNM frequency is recovered. These 
observations have important implications for interpreting the highly damped 
QNM's. 

Of course the QNM's cannot, by definition, probe the black hole interiors 
directly. However, for each independent scale associated with the black hole 
exterior, there is an appropriate range of dampings over which the QNM spectrum 
is sensitive to the structure of the exterior metric at that scale. It is perhaps worth 
elaborating this point a bit further. The equations governing the propagation of 
black hole QNM's can invariably be reduced to a scalar equation of Schrodinger 
form: 

^r + iW = o, (1.3) 
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where the "potential" R(r) depends explicitly on the form of the metric, and on 
the QNM frequency u. For example, for pure electromagnetic and gravitational 
perturbations in the background of a 4 — d R-N black hole, the potential can be 
approximated by 



- t_l/i {2GMr + l±^ 2GMq ^l sq j (1 . 4) 

when q « GM . Here j is the spin of the perturbation, which is for scalar, 1 
for electromagnetic and 2 for gravitational perturbations. In the case of infinite 
damping, i.e. 

GM\u\ » » 1 , (1.5) 

the analytic techniques are only consistent if one retains the last term in R(r), 
despite the fact that exterior to the horizon, r > 2GM, this term is much less 
than those that depend only on M. That is 

4/6 4 

q ^ < 77^1 « 1 • (1-6) 



(GM) 2 /r 4 {GM) A 

This is what is meant by the statement that for suitably large damping, the 
corresponding QNM spectrum is sensitive to the small scale structure of the 
spacetime. On the other hand, for intermediate damping, as in ( 11. ip . one can 
safely ignore all terms containing q so that the spectrum is purely Schwarzschild 
and contains no information about corrections due to the presence of electric 
charge and the inner horizon. Similar arguments can be made for the Kerr case 
as well. 

The above discussion leads one to ask about the role of the classical singularity 
in determining the highly damped QNM's. At some scale, presumably the Planck 
scale, quantum gravity corrections will become important. A great deal of work 
has recently been devoted to the question of how, or indeed if, the singularity 
at the center of black holes gets resolved by quantum gravitational effects. Most 
of this work has been in the context of loop quantum gravityjSj-[I3]. One can 
therefore ask how the resulting quantum gravity motivated corrections to black 
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hole spacetimes exterior to the horizon affect the highly damped QNM's. Or 
to put it differently, can one even in principle obtain information about the 
quantum gravity induced small scale structure of a quantum corrected (QC) black 
hole by measuring its highly damped QNM's? 

The purpose of this paper is to examine this question in the context of the 
particular QC black hole spacetime derived by Peltola and one of us in [H] using 
an effective polymer quantized theory of the homogeneous Schwarzschild interior. 
This QC black hole has several properties that make it particularly well suited for 
such an investigation. First of all, it is possible to explicitly write down a relatively 
simple analytic expression for the metric. This is a necessary condition for doing 
the analytic calculations of the highly damped QNM's. Secondly, the solution 
is completely regular everywhere, but has only a single horizon. The exterior 
is asymptotically flat and at large distances takes the usual Schwarzschild form. 
Using the language of the maximally extended Kruskal spacetime to describe 
the interior, one finds that the throat of the Einstein-Rosen wormhole contracts 
to a minimum "bounce radius" before expanding to infinity. Thus the interior 
asymptotes to a Kantowski-Sachs type homogeneous but anisotropic spacetime. 
The singularity is avoided because, in effect, r = is not in the manifold. It is 
therefore intriguing to see how the non-singular nature of the interior affects the 
highly damped QNM's as (in principle) seen in the exterior. 

It should be noted that there are two commonly used methods for the analytic 
(as opposed to numerical) calculation of asymptotic QNM's. The Motl-Neitzke 
(MN)[l5] method was first used to obtain analytic expressions for the highly 
damped QNM's of Schwarzschild black holes in four dimensions and higher. These 
results were later confirmed using a slightly different method by Andersson- Howls 
(AH) [2]. Both methods have since been applied to a large variety of black hole 
spacetimes . In some cases, such as the Gauss-Bonnet black hole, it turns out 
that the analytic structure of the solution makes it difficult to implement the AH 
procedure [TB]. This is related to the presence of a branch point singularity in 
the Schwarzschild form of the metric, which appears in the complex plane. The 
QC black hole that is examined here, despite being non-singular, also has such a 
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branch point singularity in the Schwarzschild form of the metric and a preliminary 
examination points to potential difficulties with the AH method in this case as 
well. We will show that the MN method can nonetheless be implemented in a 
rather straightforward fashion. 

The structure of the paper is as follows. In section 2 we review briefly the 
AH and MN methods for the Schwarzschild black hole, introducing the notion of 
Stokes and anti-Stokes lines for the uninitiated. Section 3 is devoted to a review of 
the QC black hole of [H], including a presentation of the Stoke/anti- Stokes lines 
for this well as a demonstration of why the AH method is problematic. 

In section 3, we also present our calculation of the asymptotic QNM's for the QC 
black hole using the MN method. Finally, we close in section 4 with conclusions 
and prospects for future work. 



2 Highly Damped Quasinormal Mode Calcula- 
tions 

2.1 General Method 



In order to illustrate the method in its simplest form we consider the QNM's 
of a massless scalar field propagating in a Schwarzschild black hole background. 
The equations for the gravitational wave modes, after the appropriate tensor 
decompositions [17], are qualitatively similar. The Klein-Gordon equation we wish 
to consider has the form 



where g^ is the metric and g is it's determinant. In a completely general, spher- 
ically symmetric, static spacetime, the line element can be written as 




(2.1) 



ds 2 = -f(r)dt 2 + gir^dr 2 + r 2 dtt 2 . 



(2.2) 



We separate variables (r, t) as follows 



^(r,t) = e 



(2.3) 
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where F = \[fg and u = u r + iui is in general a complex frequency such that 
bj r > and uji < 0. The spatial part of the wave function ip(r) obeys 

d 2 ib , . , 
+ R(r)i> = , (2.4) 

where the effective potential R(r) is 

K[r) ~ F 2 2 F AF 2 r F ■ ^ 
Note that for non-extremal horizons F has a simple zero at the horizon location. 
Moreover, for asymptotically flat spacetimes, F — > 1 as r — > oo so the potential 
R — > u 2 . The asymptotic solutions are therefore ingoing and outgoing spherical 
plane waves, as expected. 

The QNM's are defined as usual to be solutions to the wave equation that are 
purely ingoing at the (outer) horizon and outgoing at spatial infinity. In the limit 
that the magnitude of the imaginary part of u is large (this will be made more 
precise in what follows), the only terms in R(r) other than the term containing 
uj 2 , which are relevant to the calculation, are those that dominate near any poles 
that may be present in R. 

In the case of the Schwarzschild solution, the line element is given by 

ds * = _(i _ ^-)dt 2 + (1 - ^)-^r 2 + rW . (2.6) 

In this case there are poles in R(r) only near the horizon and the origin. In the 

infinite damping limit, the relevant part of the potential is simply 

u 2 4? 2 - 1 
*M « - ' (2 ' 7) 

where j = for scalar perturbations. 

The analytic calculation is based on the Wentzel-Kramers-Brillouin (WKB) 

approximation, which gives the two linearly-independent solutions to be 



±z J Q{r')dr' 



{2.1 



where t is one of the zeros of the function Q(r). Generally, one needs to take 

Q 2 (r) = /2(r)-£-A-, ^ 



4 ( r -P)2 



where P are the location of the poles of R. We define the square root of Q 2 in 
the WKB solution so that the positive exponential term (i.e. f±) is outgoing at 
infinity. 

The basic idea of both the AH and MN analytic methods is to find a closed 
contour around the pole at the event horizon in the complex r— plane along which 
the WKB approximation remains valid and evaluate the change in the phase of 
the above expressions (the monodromy) along this contour. If the contour can 
be suitably connected to r = oo, one can impose the boundary condition that 
the solution be purely outgoing (f\) there as a "starting point" and see how 
the solution changes as one goes around the contour. If in addition the chosen 
contour can be deformed without obstruction to a contour close to and encircling 
the event horizon, one can impose the second boundary condition that the solution 
is purely ingoing at the horizon to evaluate the monodromy along this contour. 
Requiring the two monodromies to be equal gives the WKB condition on the 
QNM frequency and determines the spectrum. A useful concept in choosing the 
contours is the anti-Stokes line, which is a line along which Q(r)dr is purely real 
and thus our WKB solutions are periodic functions. The anti-Stokes lines far 
from the origin (r = 0) and two suitable contours (path 1 and path 2) for our 
calculations are shown in Fig. [TJ Note that path 2 starts at point a. Then it 
moves to point c and after encircling the horizon it returns back to point a. This 
was the loop chosen by Andersson and Howl£[ The other contour, path 1, was 
the one used by Motl and Neitzke. It starts at point a in Fig. HJ Then, it moves 
to point d and follows the line shown in the figure to infinity and circles around at 
infinity in the counter-clockwise direction to meet the line in the upper half plane 
and finally it returns back to point cq. Since path 1 can be contracted to path 
2 without encountering any poles, one expects them to yield equivalent results, 
which indeed they do. 



1 Note that Andersson and Howls encircle the horizon in the clockwise direction in [2], but 
path 2 runs in the counter-clockwise direction. This difference should not and does not affect 
the results. 

2 Note that, in [15] . Motl and Neitzke start their loop at point d and rotate in the clockwise 
direction along path 1. The reason for this difference is the fact that Motl and Neitzke assume 
that perturbations depend on time as e , which means LOi > 0. In this paper we assume that 
perturbations depend on time as e _lwt , which means Wj < 0. 
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Fi gure 1: A simplified schematic illustration of anti- Stokes lines far from the origin for the 
Schwarzschild spacctime in the complex r-plane. The filled circles are the poles of the function 
Q(r) at the event horizon (r = r/j) and the origin (r = 0). The dashed lines represent the paths 
that we follow along anti-Stokcs lines to determine the QNM spectrum. 

As long as one stays on an anti-Stokes line the relative weights of the two 
WKB solutions do not change. So if the solution is purely f±, it remains j\ 
as long as one does not leave the anti-Stokes line. The problem, as illustrated 
by the two contours in Fig. [1] is that one cannot find a suitable closed contour 
that remains on a single anti-Stokes line, since in general they emanate from the 
neighborhood of poles in R(r). The trick is to know how the solution changes as 
one "jumps" from one anti-Stokes line to another in the neighborhood of a pole. 
Motl and Neitzke[15] solved the approximate wave action analytically near the 
poles in order to calculate the monodromy along the contours in their vicinity. 
Andersson and Howls [2] clarified the MN calculation by using the so-called Stokes 
phenomena to calculate the monodromies in the vicinity of the poles. In all 
previous calculations in the literature, both methods gave identical results, as 
they must. 

We will now describe both in the context of the Schwarzschild black hole 
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before proceeding to apply them to the non-singular black hole of reference [2]. 



2.2 The Andersson-Howls Method Applied to Schwarzschild 
Black Holes 

The AH method requires a detailed knowledge of the analytic structure of the 
function Q(r) near the pole at r = 0. For this, one needs to determine the location 
of the zeros of Q in the complex r-plane. It is important to do this as it is from 
the zeros that the Stokes and anti-Stokes lines emanate. There are three of each 
emerging from each zero in the complex plane, being initially spaced 2-7r/3 radians 
apart from one another. We also need to introduce a branch cut from each zero 
to ensure that Q remains single-valued. These cuts need to be introduced in a 
way that they do not affect our analysis. Specifically, they should not intersect 
path 1 and path 2 in Fig. [2J The Stokes lines are the lines in the complex plane 
where Q(r)dr is purely imaginary (and thus one of the WKB solutions \ blows 
up exponentially as one moves away from the zero of the function Q(r) at t 
along Stokes lines). When one moves from one anti-Stokes line to an adjacent 
anti-Stokes line in the vicinity of a zero, a Stokes line has to be crossed. In this 
case, Stokes phenomenon dictates that the WKB solution, which is dominant on 
the Stokes line, be replaced by the dominant WKB solution ±i times the sub 
dominant term. The positive sign is used when one crosses the Stokes line in the 
anti-clockwise direction and the negative is used for the clockwise crossing. Note 
that such simple change, i.e. the addition of ±i times the sub dominant WKB 
solution to the dominant one, only holds when the Stokes line emanates from 
the zero, which is used as the lower limit in the phase-integral that appears in 
the WKB solutions ff\. This means that in order to construct an approximate 
solution valid in different regions of the complex plane we need to switch the 
lower limit of our phase integral between zeros. This requires the evaluation of 
integrals of the type 



where tj and tj are two zeros of Q. To derive a WKB condition on QNM's, we 
can take either path 1 or path 2. In both paths, we start on line labeled a on 




(2.10) 
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Figure 2: A complete schematic illustration of Stokes (dashed) and anti-Stokes (solid) lines 
for the Schwarzschild spacctimc in the complex r-planc. The hollow circles are the zeros and 
the filled circles are the poles of the function Q(r). The poles are located at the event horizon 
(r = Th) and the origin (r = 0). The paths {path 1 and path 2) that we can follow along 
anti-Stokes lines to determine the QNM spectrum are represented by thin dashed lines with 
lighter colors. 

which we know what the solution is due to the boundary condition at infinity: 

*a = /f l} , (2.11) 

where the superscript t\ indicates that the WKB phase integral is to be evaluated 
from the zero t\ shown in Fig. [2j In the case of path 1, we can move along anti- 
Stokes lines, while applying Stokes phenomenon when we jump from one anti- 
Stokes line to the other near zeros, to the line labeled / that extends to infinity 
in the lower half plane. On line /, we get the solution 

V f = -(e" 3i7 + e" 47 + e i7 )/f 4) + i{e~^ + e" 4 ' 7 + e 4 ' 7 + e 3i7 )/ 2 (il) , (2.12) 

where 

7 = 731 = 732 = 742 , (2.13) 
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which can be evaluated analytically in the limit r — > 0. In this limit, we can use 
Eqs. (J22D, O, and fl2~T0|) to show that 



731 



t 3 



r 2 u 2 j 2 ~ 



2 dr = |tt = 7 . (2.14) 



-2GM) 2 r 4 

We can now close path 1 at infinity, where f\ is dominant. After returning 
back to line labeled a, we get 

ty a = -(e" 3i7 + e~ i7 + e i7 )e i(rccw+7) /i (tl) + sub dominant term , (2.15) 

where ccw stands for counter-clockwise and 

Fcoa = <P Qdr = 2m Res (Q) = AniuGM. (2.16) 

J r=2M 
ccw 

We now can use the fact that 

# s = e-^tta , (2.17) 

where e~ lTccw is the monodromy of our loop according to the boundary condition 
at the event horizon to get the WKB condition 

e -2iv ccw = e s™GM = _ x _ 2cos (2 7 ) = -1 - 2cos(j7r) (2.18) 

on the QNM frequency u. The disadvantage of path 1 is that, as pointed out by 
Motl and Neitzke, it can only be trusted near infinity for the dominant mode, so 
that the sub dominant mode must be dropped. While this is sufficient to give the 
correct condition on the frequency, one gives up the ability to do a consistency 
check using the monodromy of the sub dominant mode. The AH contour [path 
2), however, does not suffer from such disadvantage. In the case of path 2, after 
circling around the horizon, we get back to the line labeled a, where we get 

= _(1 + e 2i7 + e- 2i ^)e iT f[ tl] - i(l + e 2 * 7 ) [e~ lV + e ir (l + e 2 * 7 + e" 2 * 7 )] / 2 (4l) . 

(2.19) 

Comparing and \I/ a tells us that the coefficient of the sub dominant WKB 
solution f 2 has to be zero. This gives us the same WKB condition as in Eq. 
(12.181) . We now can use (12.181) to show that the coefficient of the dominant WKB 
solution fx is equal to the correct monodromy presented in Eq. (I2.17p . This is 
what we mean by consistency check. 
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2.3 The Motl-Neitzke Method Applied to Schwarzschild 
Black Holes 



An alternative means of determining the QNM spectrum is given by Motl and 
Neitzke. For this method the simplified version of the anti-Stokes line topology 
in the infinite damping limit, shown in Fig. (TJ is sufficient. In other words we 
can ignore the "fine structure" of zeros surrounding the pole at the origin, as well 
as the details of the Stokes lines. The basic idea is to approach each singularity 
in the complex plane along an anti-Stokes line closely enough so that a suitable 
approximation to the differential equation (12. 4p can be solved analytically. In 
all cases that have been considered so far, the solution near the poles is well 
approximated by a Bessel function. This approximate analytic solution is then 
used instead of the Stokes phenomenon to "jump" between anti-Stokes lines near 
a pole. 

Our task then is, essentially, to "walk" our solution around a closed contour 
such as path 1 or path fusing the WKB approximation along the anti-Stokes lines 
and the approximate analytic (Bessel function) solution near the poles. 

We now calculate the monodromy around path 1 using the MN method. Note 
that we have chosen the branch cut from the origin to point along the negative 
real axis. We start at the point a in Fig. [1], such that r a << 2GM which is near 
enough to the pole at r = so that Eq. (12.41) can be well approximated by 



d 2 m 



dr 2 



+ 



r 2 u 2 4j 2 - 1 



-2GM) 2 4r 4 



* = . (2.20) 



The general solution to the above is: 

* = dy/wrJi (z) + C 2 V^J-i (z) , (2.21) 

2 2 

where 

9 

r uj 

Z = —7777 2.22 

and J± v are Bessel functions of the first kind. We now need to relate this general 
solution, near the pole at r = along the anti-Stokes line labeled a to the 
outgoing and ingoing WKB wave solutions ( 12. 81) . Sufficiently far from the pole at 

(0) _ I-2GM n iz n „A f(°) ~ / —2GM 



r = these are well approximated by f{ ~ y ^ etz an d / 2 ~ \l ^ e " 
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respectively. Next we evaluate the asymptotic form of the Bessel functions. To 
do this, we use the fact that 

J±„(z) = z ±v <j>{z) , (2.23) 

where <fi is an even holomorphic function. This means we can write 

J±„ (e«e-*z) = e ±w6 J± v {e~ ie z) , (2.24) 

where 9 = argz. Note that e~ l6 z = \z\ is positive and real. As long as we stay 
away from r = at a distance where 



AGM 

|r|»W^— , (2.25) 

we can assume that \z\ » 1. This allows us to use the asymptotic behavior of 
the Bessel function 

/ 2 / 7T 7T \ 

J± v {x) — y \j — cos [x =f v- — — ) for x » 1 (2.26) 
V 7TX V 2 4/ 

to show that 

s/^J± v {z) = ^e ±lu9 J ±l/ (e- l0 z) 

M -e^ a± cos(e-^ - a±) , (2.27) 



where a± = ^ ± z/|. If we choose the branch — 7r < < 7r, the argument of z 
along the anti-Stokes line labeled a is 6 = 2ir. After substituting 2ir for 9 in the 
above equation, we find the wave function along anti-Stokes line a in terms of the 
ingoing and outgoing WKB modes, which is 

m a = (C ie 3ia + + C 2 e 3ia -) fi°\z a ) + (C ie 5 ^ + C 2 e 5ia -) f { 2 °\z a ) . (2.28) 

The boundary condition at infinity tells us that, along anti-Stokes lines that 
extend to infinity, the solution must be dominated by the outgoing mode. In 
other words 

de 5ia+ + C 2 e 5ia - = . (2.29) 
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We now repeat the above expansion near the point b. The anti-Stokes line b is 
rotated relative to line a one by —ir/2 in the r-plane, which means z b = e~ m z a . 
Therefore, 9 should be replaced by n in Eq. f)2.27p . This yields, after a few steps: 

v& 6 = (C ie 3 ^ + C 2 e 3 ^) ff\z b ) + (C ie tQ + + C 2 e ia ~) ff\z b ) . (2.30) 

Repeating this to get to c and then finally to d, we find that ultimately the 
— 3ir/2 rotation in the complex r-plane that takes us from the point a on the first 
anti-Stokes line to the point d on the last anti-Stokes line yields 

V d = (C ie - 4Q + + C 2 e~ ia -) f[°\z d ) + (C ie - 3ia + + C 2 e~ 3ia -) f ( 2 °\z d ) . (2.31) 

The trick now is to realize that we can extend the expression for ip d and evaluate 
it in terms of z a by using the following integral expressions that are understood 
to be evaluated along the large loop at infinity: 

r Vd r r a 

e ±iZd = exp(±z / Qdr) = exp(±z ® Qdr) exp(±z / Qdr) = e ±ir ^ e ±iz * 
Jo J Jo 

ccw 

(2.32) 

so that 

V d = (C ie - 4Q + + C 2 e~ ia -) e ir ^/ 1 (0) (^)+(C 1 e- 3iQ! + + C 2 e~' Ma -) e~ iV - ff> \z a ) = V- a . 

(2.33) 

These two expressions ^> a and \l/ a at z a must differ by the appropriate monodromy 
of the WKB phase, which can easily be evaluated by invoking the boundary 
conditions. In particular, the path 1 can be shrunk without obstruction to a 
small loop encircling the horizon. Given that the solution near the horizon is 
ingoing, it must be true that 

m a = e~ iFccw ^ a . (2.34) 

The above condition, applied to the coefficient of the dominant solution, yields 
the following equation: 

(de- ia+ + C 2 e~ ia -) = e~ iVc ™ (Cie 3ia+ + C 2 e 3ia ~) . (2.35) 
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Equations (12.291) and (I2.35j) are two linear equations on the two initial coefficients 
Ci s 2, which can only be satisfied if the corresponding determinant vanishes. This 
yields precisely the WKB condition ("12 . 1 8[) as before. 

Note that with the large loop (path 1) we cannot look at the condition on 
the sub-dominant mode, so it does not yield any further conditions. However, 
by looking at the same calculation using path 2, one can find two different sets 
of linear equations using the dominant and sub-dominant modes. Both of these 
sets of equations in turn should give the same WKB condition in f 1 2 . 1 8 j) . This 
provides a very useful consistency check of the calculation, as will be seen in more 
details when we consider the more complicated case of the QC black hole. 

3 The Quantum Corrected Black Hole 

The QC black hole we wish to consider was derived in [H] using a loop quantum 
gravity inspired quantization scheme similar to those that were applied to black 
hole interiors by a variety of authors[6l [101 HH [12]. The basic idea is to start with 
homogeneous slicings of spherically symmetric spacetimes in Einstein gravity. 
Polymer quantization 0[8] is then applied to the resulting Hamiltonian to produce 
a set of effective, QC equations for the geometrical variables. These equations 
are solved to find the QC metric in the Schwarzschild interior. In contrast to the 
previous work, reference [H] only applied polymer quantization to the phase space 
variable corresponding to area. This resulted in a simple analytic expression for 
the interior metric that could be analytically continued to the exterior in order 
to construct the complete QC spacetime. The metric for the QC black hole in 
Schwarzschild-like coordinates is: 

^Ur + r W , 

(3.1) 

where k is a new (quantum) length scale, and e can be 1 or -1. 

The above metric has intriguing properties. It contains a single horizon at 
r/i = a/ (2GM) 2 + k 2 so that the analytically continued spatial slices are topo- 



, o / r ¥ 2GM\ , 2 
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logically equivalent to the Einstein-Rosen wormhole. As with Einstein gravity, 
the radius of the throat of the wormhole is time dependent. The maximum is r^, 
but instead of shrinking to a singularity in the black hole interior in finite time, 
the throat radius reaches a minimum at r = k, before re-expanding indefinitely 
to a Kantowski-Sachs spacetime. The resulting complete, single-horizon, non- 
singular black hole spacetime scenario is therefore reminiscent of past proposals 
for universe creation via quantum effects inside a black hole [9]. 

To find the QNM's we insert the QC metric ((HUD into Eq. (EH)- Using the 



substitution x — Vr 2 — k 2 and extending the domain of x to (—00 < X < 00) 
(which effectively negates the need for e), we obtain the differential equation 

+ R(xW = , (3.2) 

where ip is our rescaled wavefunction and our new potential function R(x) is given 
by 

~ u 2 (x 2 + k 2 ) k*-6 X 2 k 2 , GM 



( X - 2GM) 2 A X 2 (X 2 + k 2 ) 2 (x 2 + k 2 )( X ~ 2GM) 
GM x {GM-x? 



X 2 (X-2GM) ( x 2 + k 2 )(x-2GM) ' X 2 (x ~ 2GM? ' (3 ' 3) 
We immediately notice that this function has poles at x = and x = ±ik. It also 
has a pole at the point x = 2GM, which corresponds to the event horizon. Note 
that in the infinite damping limit, the relevant part of this potential is simply 

- _ u\ X 2 + k 2 ) 1 7k 2 
{X) ( X - 2GM) 2 x 2 4(x 2 + k 2 ) 2 ' 1 ' 

This is because the last two terms on the right hand side of the above equation 
dominate all the other terms near the poles whereas the term containing u cannot 
be neglected elsewhere. The two linearly-independent WKB solutions to this 
equation are given by (12. 8p with 

Q 2 (x) = R{x) - —, — ^-7-2 -j~2- 77 — —772 - 77 7777772 • ( 3 - 5 ) 

The last term above can be neglected for large damping rates. 

Finally, we want to point out that in an intermediate damping region of 

1 « GMU « (3.6) 
k z 
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Figure 3: A complete schematic illustration of Stokes (dashed) and anti-Stokes (solid) lines for 

the QC spacetime in the complex x-plane when |cj| < (left) and |w| > ^ 2G ^P (right). 

The hollow circles are the zeros and the filled circles are the poles of the function Q(x)- 

the potential ( 13 .4[) reduces to the Schwarzschild potential given in Eq. (12. 7p . 
Therefore, in the above intermediate region of the spectrum, the Schwarzschild 
value of ln(3) is expected to appear in the real part of the QNM frequency u. 

3.1 The Andersson-Howls Method Applied to the Quan- 
tum Corrected Black Hole 

In Fig. |3l we show schematically the topology of Stokes and anti-Stokes lines in 
two different damping regions for the QC black hole under consideration. The 
topology on the right appears in higher damping rates compared to the topology 
on the left. Fig. [3] is plotted based on numerically generated figures of Stokes and 
anti-Stokes lines, some of which are shown in Figs. H] and [5j After investigating 
the numerically generated topology of Stokes and anti-Stokes lines for different 
values of M and k, we conclude that in the topology shown on the right in Fig. [3] 
the anti-Stokes line that connects the poles at \ = ±ik misses the pole at \ = 
by a minimum distance of ~ -^~ r When we lower the damping, the zeros of 
the function Q(x) shown with hollow circles in Fig. [3] move away from the poles. 
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Figure 4: Numerically generated anti-Stokes lines for the QC black hole spacetime in \~ 
coordinate, showing the change in topology as one goes from \u>\ < '" 2G ^P to \ui\ > G ffl ■ 
The figure on the left shows the structure of the lines near the origin of the complex plane 
for k = 0.1, 2GM = 10 and ui = -ilO 5 . In the figure on the right, k = 0.1, 2GM = 1 and 
u> = — z5000. The filled circles are the poles of Q located at ±ik and and the smaller hollow 
circles are the zeros of Q. 

Using Eqs. (13. 4p and (13. 5p . it is easy to show that the distance between the zeros 
and the pole at x = is ~ Equating this with the minimum distance that 

the anti-Stokes line misses the pole at \ = leads us to an approximate location 
of the damping rate where the topology changes from the one on the left to the 
one on the right in Fig. |3j This damping rate is 

k| « M w fc3 ■ (3-7) 

As we have seen during the Schwarzschild case, one ideally needs a closed 
anti-Stokes contour around the pole at the event horizon, which connects to an 
outgoing anti-Stokes line. Critically, however, we require that such a contour 
never intersects with any pole. A cursory analysis of Fig. |3] shows that no such 
contour exists in either topology. Therefore, this method appears to be problem- 
atic for this case. 
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Figure 5: Numerically generated anti-Stokes lines for the QC black hole spacetime in \~ 
coordinate for k = 0.1, 2GM = 10 and u> = — HQ 9 . The filled circles are the poles located at 
±ik and 0. 

3.2 The Motl-Neitzke Method Applied to the Quantum 
Corrected Black Hole 

In this section, we evaluate the highly damped QNM's of the QC black hole 
under consideration using the MN method. In Fig. [61 we show the topology 
of anti-Stokes lines far from the poles, which is sufficient for applying the MN 
method. As before we construct a simplified contour (path 1 or path 2 in Fig. [6]), 
where our task is to "walk" our solution from an anti-Stokes line, which extends 
to infinity, around the pole at the event horizon and back to the same anti-Stokes 
line without crossing any of the branch cuts. As it was pointed out in the last 
section, at some damping rate above ^ 2G ffi , the anti-Stokes line that connects 
the poles at \ = ±ik will miss the pole at \ = 0. However, for physical reasons, 
we are interested in the limit where k « GM. In this limit, the function Q(x) 
can be approximated in the vicinity of the poles at \ = and x — ±ik, i.e. when 
Ixl ~ k, by 

The above function leads to the topology shown in Fig. [61 where the anti-Stokes 
line that connects the poles at \ = ±ik goes through the pole at x = 0. This 
is evident in Fig. [5j where we have numerically generated the topology of anti- 
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Figure 6: A simplified illustration of anti-Stokes lines far from the poles for the QC black hole 
when k « GM . The outermost pole represents the black hole's event horizon, with the three 
poles at left corresponding to those at ±ik and 0. The dashed lines are the suitable paths that 
we follow to determine the QNM spectrum. 

Stokes lines using the function Q(x) in (13.51) . in which the approximate R(x) m 
(13.41) is used, for the case where k << GM (k = GM/ 50). In the rest of this 
paper, we do the calculations assuming that k « GM. The procedure to solve 
for the QNM frequencies, for the topology shown in Fig. |6l is as follows: 

We evaluate our solution in the immediate vicinity of the pole at x — ik. In 
this region, differential equation (13. 2p can be expressed as 



d 2 j) 
dx 2 



+ 



u 2 (2ik)(x-ik) 7 



-2GM) 2 16( X -ik) 2 



ij> = . (3.9) 



This is a Bessel differential equation with a solution of the form 



An An 
^ = A + ^ — (x- ik)J v (Z) + A.\j — (x - ik)J_ v {Z) , (3.10) 

where A± are constants, v = y -jj, J±v are Bessel functions of the first kind and 



2 / W2^\ / 

3 ^GAf J y 
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Z is multi-valued. To make Z single-valued, we choose a branch cut in which 
Z a = \Z\e 2l7T along the anti-Stokes line labeled a. Note that the anti-Stokes 
lines labeled b and c are rotated relative to a by zl j L and =1 ^ L (in the x - pl ane ) 
respectively. Therefore, we get Zb = \Z\e l7T along the anti-Stokes line labeled b 
and Z c = \Z\ along the anti-Stokes line labeled c in Fig. [6j We now can take the 
same steps as in Eqs. (I2.23P through (12.261) to show that 



^ « A + 2e^ , I cos (e~ ie Z - a+ ) 

' cu\/2iky/x — ik 




+ A_2e'v"- / J GlK J__ cos (e~ i0 Z - a J) , (3.12) 

where 6 = argZ and a± = j (1 ± y— §Y Note that the above approximation 
of Bessel functions is only valid if \Z\ » 1. For \Z\ to be large, we need to stay 
away from the pole at x — ik at a distance where 

/ GM V 

|x-iA;| » i -j =-|w|J . (3.13) 

In addition, for our calculations to be valid, the above distance should be much 
smaller than the polymerization length scale k, which determines the distance 
between the poles at ±ik and 0. This way our calculations near the pole at ik 
will not get affected by the presence of the pole at % = 0. This leads to a lower 
bound on the size of \oj\, i.e. 

GM 

M » -gr ■ ( 3 - 14 ) 

For e~ ie = +1, we get 

+ [A + e ta ^ l+2 ^ + A_e ia -^) /f ] , (3.15) 
and for e~ 10 = —1, we get 

f/, = (A + e^( 1+ f) + A_e--( 1+ v)) / f) 

+ (A+e-^+^-v) + ^.e^-^-v)) /f > , (3.16) 
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where 



and 



/i (ifc) = J—^¥=* iZ ( 3 - 17 ) 



ft ] = J ~l G ^ =e ^ (3.18) 
are the approximate WKB solutions to the wave equation (13. 9p at a distance 

2 

>> \ 1 from the pole at ik. Using Eqs. (13 . 151) and (13 . 161) we can determine 

the wave solutions along the anti-Stokes lines labeled a, b, and c: 



4>a 


= (± 


_ e 3ia + + A 


e 3 * a ~)ff\z a )- 


H (A + e 5ia+ 


±A_e 5ia -)ft\z a 




= (An 


_ e 3ia+ + A _ 


e^-)ff k \z b )^ 


- (A + e ia+ -} 


-A^-)ft\z b ) , 


'W c 


= 


.e~ ia+ + A 


_e- ia ~) ff ] (Z c ) 


+ (A + e ia+ 


+ A_e* a -)ft ] (Z c ) 



(3.19) 

We now move to the pole at \ = 0, where the QNM wave equation takes the 
form 



with the Bessel solution 



0J 2 k 2 1 

$ = , (3.20) 



-2GMf X 



i> = B +V ^J^ (Z ) + £_ v^J_^3 (Z ) , (3.21) 

2 2 

where 5± are constants and 

cok 

To make Z single- valued, we choose the branch cut in which = \Z \e l7T along 
the anti-Stokes line labeled d. A —it rotation in the %-plane will take us to the line 
labeled e on which Z 0e = Z^^er™ = \Z \. Once again, we can use the asymptotic 
behavior of Bessel functions (I2.26P in the large Zq limit to show that 



^ = B + 2e i ^+\l-^coB(e-» Z -l3 + ) 

ujk ' 



+ B_2e i ^-J^icos(e- ie °Z -P_) , (3.23) 



22 



where 6q = aigZ and f3± = f (l ± \f— 3j. From Eq. (I3.22p . it is clear that 
the condition, which allows Zq to be large without getting too close to the other 
poles, is precisely that given in Eq. (I3.14p . This is required for the validity of the 
calculation. 

For e~ ld ° = +1, we get 

+ {B + e^ 1 ^ + B_e*-< 1+3 ^) /f , (3.24) 

and for e~ 10 = —1, we get 

V = (B + e^ 1+ ^ + B.e l ^ l+ ^) 

+ (B+e-^t 1 -^) + B-e-V-V-^) /f , (3.25) 

where 



/, (0> = l/=|jr«* (3-26) 



and 



are the approximate WKB solutions of the wave equation (13.201) when \x\ » fj^y- 
Using these equations we can determine the wave solutions along the anti-Stokes 
lines labeled d and e to be: 

1>* = + B^~) fi°\z od ) + (B + e^ + B_e#-) ff\z od ) , 

^ e = {B + e-^ + B.e^-) f[ 0) (Z 0e ) + (S +e ^ + B^~) f?\z 0e ) . 

(3.28) 

We finally move close to x — —ik, where the differential equation (13. 2p can 
be expressed as 



+ 



dx 2 

with the solution 



u 2 {-2ik){ X + ik) + 



2GM) 2 16(x + iky 



ip = (3.29) 



^ = C + y/^( X + ik)J v (Z) + C-^( X + ik)J- v (Z) , (3.30) 
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where C± are constants and 



2 / uiy/— 2ik \ 

Z=o Ms™ )(x + ik) 3/2 . (3.31) 



3 V -2GM 

To make Z single- valued, we choose a branch cut in which Zf = \Z\e lT: along 
the anti-Stokes line labeled /, Z g = \Z\ along the anti-Stokes line labeled g and 
Zh = \Z\e~™ along the anti-Stokes line labeled h. Following similar steps as in 
previous cases, we find 



V>/ = 




,e 3ia+ + C_ 


<■'"•■ )/; • 


■ k \z f ) 


+ (C + e ia+ +C_e m -) ft 




Ipg = 




_ e - ia+ + c 




- lk \z g ) 


+ (C + e ia+ +C.e ia -) / 2 ( " 


- ife) (z 9 ) , 


Iph = 




_ e - ia+ + c 


-e- ia -)f { { 


- lk \z h ) 


+ (C + e- 3ia+ + C_e~ 3ia - 


)ft ik \z h 



(3.32) 



where 



and 



= J—^¥=e a (3.33) 



/<-*> = _ (3.34) 

are the approximate WKB solutions to the wave equation (I3.29P at a distance 

2 

>> ^ ^|^| J 3 from the pole at —ik. 

The solutions do not change in character along anti-Stokes lines. This means 
that, since the anti-Stokes line labeled a extends to infinity, we can impose the 
boundary condition at infinity on ip a . This gives us the condition 

A+e 5ia+ + A_ e 5ia - = . (3.35) 

Also, ip c and ^> which are on the same anti-Stokes line, have to be equal. Since 
we can write 



^ = [A + e~ ia+ + A_e~ ia ~) e^ d ff\z d ) + (A + e ia+ + A_e ia -) e~ 5 ° d / 2 (0) (Z, 



d 



(3.36) 
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where 



Cd L Q X L { X- 2GM ) dX ~ L { -2GM ) dX 8GM 

(3.37) 

appears as a result of changing the lower limit of the phase integral of the WKB 
solutions from ik to 0. We can now compare ( 13.361) with ipd given in Eq. ( I3.28p . 
which gives the following two conditions: 

(A + e~ ia+ + A_e~ ia -) e^ d = B + e 3ip+ + B_e 3ifi - , (3.38) 

(A + e ia+ + A_e ia ~) e~^ d = B + e ip+ + B_e ip ~ . (3.39) 
The same thing can be done for ip e and ijjf, which results in two more conditions: 

{B + e- ll3+ + B.e-tf-) e i5 ^ = C + e 3ia+ + C_e 3ia - , (3.40) 

(B+e if)+ + B_e lfi -) e~ iS ^ = C + e ia+ + C.e ia ~ , (3.41) 

where 

5 ef = / Qd X = 5 cd . (3.42) 

So far, all the conditions above apply to both path 1 and path 2 in Fig. ((6]). In 
the case of path 1, since the anti-Stokes line h extends to infinity, we can impose 
the boundary condition at infinity to iph, which gives us the condition: 

C + e~ 3ia+ + C.e~ 3ia - = . (3.43) 

This condition however does not impose any WKB condition on the QNM fre- 
quency. Another condition, that we can get, is by moving from point h to infinity 
along the anti-Stokes line that extends to infinity in the lower half of the complex 
X-plane. Then we can rotate at infinity, where we know the solution due to the 
boundary condition, to the other anti-Stokes line which extends to infinity in the 
upper half of the complex plane. Finally we can move along this anti-Stokes line 
all the way back to point a, where we get 

ijj h = (C + e~ ia+ + C-e~ ia -) / 5h «fi k \z a ) + sub-dominant term = • (3.44) 
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For this path the coefficient of the sub-dominant solution f2 is not reliable. Here 

rik 

ha= / Qd X (3.45) 



J —ik 

along the path to the right of the event horizon. By returning back to point a, 
we have completed a loop around the event horizon. Therefore, in matching ijj a 
and ipa we need to account for the monodromy of path 1. In other words 

Va = e~ iA -^ a , (3.46) 

where A ccw is the integral of Q along a contour encircling the pole at x — 2GM 
in the counter-clockwise (ccw) direction: 



A ccw = I Qd X « 2m Res \ ^ X 1^^ \ = 2vr W4G 2 M 2 + k 2 
J X =2GM \ x ~ 2GM } 



(3.47) 



Equation (I3.46P gives us the final condition 

(C + e~ ia+ + C_e~ ia -) e i^ e -^c d +s ef ) = (A +e 3ia + + A_e 3ia ~) e ~ iA ^ , (3.48) 

where we have replaced Sh a with A ccw — 5 c d — 5 e f. 

In the case of path 2, we first need to match the solution at g with the solu- 
tion at b, which are on the same anti-Stokes line. This matching results in two 
conditions: 

(C + e~ ia+ + C^e~ ia ~) e i*cc W e -i(5 cd +s ef ) = A +e 3ia + + A_e 3ia ~ , (3.49) 
(C + e ia+ + C-e ia ~) e - iA ^ e KSc d +s e} ) = £ +e ia + + A_e ia ~ . (3.50) 
Here we have used the fact that 

/ik 
Qdx = A cc „, - 5 cd - 5 e f , (3.51) 
■ik 

where 5 g b is taken along the anti-Stokes line to the right of the pole at the event 
horizon. We finally can return to anti-Stokes line a where we can impose the 
boundary condition at infinity to the coefficient of the sub dominant solution to 
get 

A+e 5ia+ + A_e 5ia ~ = . (3.52) 
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We can also impose the monodromy condition that "0a 



-iA c 



' J ip a on the coeffi- 



cient of the dominant solution to find 



A + e 3ia + + A_e 3ia - = (A + e 3ia+ + A_e 3ia -) e~ iL 



(3.53) 



Finally, we can use Eqs. (137351) . (IH755|> . ( 1373T)j) . ( 1374Dj) and (EH]) together with 
either (I3.48P from path 1 or with the combination of Eqs. (I3.49p . (I3.50p and (I3.52p 
or (I3.53P from path 2 to get a WKB condition on u. All the above combinations 
should result in the same WKB condition on the QNM frequency, which they 
do. This provides a valuable check on the validity of our results. For the case of 
path 1, in order to have a nontrivial solution to Eqs. (I3.35p . (I3.38p . (13. 39 p . ( I3.4(jp . 
(I3.4ip and (I3.48p . the determinant of the matrix 



,5i Q4 



n 5i a- 



e ia +e -i5 e ia- e ~i8 






























3i a+ 
D i a_)_ 



_g— ja+g2iA CCUI( 



— e 



-2i<5 



has to vanish. This leads to the WKB condition 







3i ct- 



— e 



(3 



e -2iA ccw + Kie -*iS + K2e -2iS + K 3 = 



(3.55) 



where 



5 
K 3 



4 cos 2 (a + — a_) , 

4cos(a + — a J) cos(/3 + — /?_) 

4cos 2 (a+ — a_) — 1 . 



(3.56) 



Since A ccw m 4-kiuGM + 45, the above equation can be written as 



TTUjk^ ZlTUjkP' TTUjk^ 



(3.57) 
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We now take y = 8iruGM to get 

e y + K ie - 29y + K 2 e~ 39y + K 3 e' 4gy = (3.58) 
where g = 32 q2 M 2 ■ We then take y = ( — ir), which leads to 

e c cos 77 + K ie ~ 29C cos{2gr]) + K 2 e~ 39<: cos(3#r/) + K 3 e~ 4aC cos(4#r?) = (3.59) 
and 

- e c sin r/ + K x e~ 2 ^ sin(2^r/) + K 2 e~ 3gC sin(3^r/) + K 3 e~ i9C sm(Agr]) = .(3.60) 

From the last two equations, we see that for the solution to be periodic in 77, we 
must require simultaneously 

rj — > t] + 2n7r and g(rj + 2nir) = grj + 2mir , (3.61) 

where n and m are integers. In other words, periodicity in the damping term is 
only possible if 

m 

9 = - ■ (3.62) 
n 

Note that since g = 32 g^m 2 ' we rec L mre — 1- 

For m and n integers we can introduce z = e y ^ n to obtain 

z n+4m + RiZ 2m + g^m + ^ = Q , (3.63) 

which has n + 4m roots. (Note that in z n+im above, one can neglect 4m since 
n » m. This is equivalent to neglecting the term AS in A ccw .) These roots lead 
to 

8ttojGM = y = n In \z\ + in arg(z) — 2npiri , (3.64) 

where p is a large integer. Using the above equation, we have plotted the QNM 
frequency spectrum in Fig. [7] for m = 1 and n = 10, 100, 1000. 
In the limit k — > 0, Eq. (I3.63P leads us to 

e toGM = -(ifi + K 2 + ITa) = -59.5392... . (3.65) 
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Figure 7: Real part of 8ttujGM versus the imaginary part for m = 1 and n = 10, 100, 1000 
from left to right respectively. 

This gives a QNM frequency spectrum of the form 

87rwC?Af = ln(59.5392...)- 27ri U + ij , (3.66) 

where I is a large integer. This is different than the Schwarzschild result due to 
the presence of the quantum length scale k in the spacetime metric. 

We, finally, would like to point out that in the intermediate damping region 
of (13. 6p . our calculations in the complex plane are done at a distance of |%| >> 
k from the origin of the complex plane. This is the reason why we get the 
Schwarzschild value of ln(3) for the real part of u. In this intermediate damping 
region, our path will not be affected by the "fine structure" of the poles near 
X = because, for example in path 1, the jump from anti-Stokes line a to line h 
in Fig. [6] is done at a distance scale much larger than k. 

4 Conclusions 

We have presented an analytic calculation of the highly damped QNM's of a non- 
singular QC black hole using the methods of Motl and Neitzke. As expected, the 
spectrum of highly damped QNM's are sensitive to the metric structure at the new 
scale, k, which determines the onset of the quantum corrections. In analogy with 
what happens in the R-N case, the real part of the highly damped QNM frequency 
does not reduce to the Schwarzschild value even in the limit that A; — ^ 0. As in the 
R-N case, this is explained by the fact that the limit \u\ — > oo does not commute 
with the limit k — > 0. The present calculation suggests that other non-singular 
QC metrics with different small scale behavior will yield a different answer in this 
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limit. This is currently under investigation. The implication is that the highly 
damped QNM's in principle, if not in practice, contain information about the 
underlying theory that is used to resolve the singularity. 

Our analysis is also interesting for the following reason. We found that, in 
analogy with the Gauss-Bonnet black hole [16], the AH method was difficult to 
implement with the given topology of Stokes and anti-Stokes lines of the QC black 
hole. This is related to the fact that both metrics have a branch point singularity. 
While these difficulties with the AH method are likely surmountable, we have 
shown that the MN method does in fact appear to produce an unambiguous 
answer in the present case in a rather straightforward fashion. This suggests 
that it may be possible to get a consistent solution for the more complicated 
Gauss-Bonnet black hole using the MN method. This calculation is also under 
investigation. 

Finally, we stress again the fact that the analytically calculated high over- 
tone QNM's probe the structure of the black hole exterior down to length scales 
determined by the inverse of the magnitude of the QNM frequency. Whatever 
corrections to the small scale structure result from the ultimate quantum gravity 
theory, QNM's which are large relative to the inverse horizon length but small 
compared to other inverse length squares will be insensitive to these corrections 
and can in principle provide information only about the horizon structure. 
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